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Abstract. We study a gas of hard rods on a ring, driven by an external thermostat, 
with either elastic or inelastic collisions, which exhibits sub-diffusive behavior {x^) ~ 
We show the validity of the usual Fluctuation-Dissipation (FD) relation, i.e. the 
proportionality between the response function and the correlation function, when the 
gas is elastic or diluted. On the contrary, in strongly inelastic or dense cases, when the 
tracer velocity is no more independent of the other degrees of freedom, the Einstein 
formula fails and must be replaced by a more general FD relation. 
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Introduction The typical scenario in diffusive problem is the so called standard 
diffusion, which is qualitatively similar to the usual behaviour in the Brownian motion, 
i.e. at large time one has 

{x{tf) ^ 2Dt , (1) 

where D is the diffusion coefficient which is linked to the velocity correlation function 
via the Kubo formula 

poo 

D= / {v{t)v{Q))dt. 
Jo 

Of course the above scenario holds if J^{v{t)v{0))dt is finite and non-zero. 

On the other hand, it is well known that, beyond the standard diffusion, one can 
have anomalous diffusion [Il[2], i.e. 

{x\t)) ~ t^" with u ^ 1/2, (2) 

formally this corresponds to have D = ooifz/>l/2 (superdiffusion) and D = ii 
u < 1/2 (subdiffusion). 

From the well estabished linear response theory, it is known that, when {x{t)) = 
in the unperturbed system, ([1]) implies a hnear drift 



x{t) ~ t , (3) 

if a small external force is applied [21 SI- In the following we will indicate with (■) the 
average in the unperturbed system, i.e. weighting states according to the stationary 
phase-space distribution and with (■) the time dependent average in the dynamical 
ensemble generated by the external perturbation. One can wonder how Eq. ([3]) changes 
in presence of anomalous diffusion, i.e. if, instead of ([I]), eq. ([2]) holds. 

The "usual" fluctuation-dissipation relation relates the mean response R{t) = 
at time t of the velocity after an impulsive infinitesimal perturbation Sv{0), applied at 
time t = 0, to the velocity autocorrelation C„(ti — ^2) = {v{ti)v(t2)): 

R{t)=C,{t)ICM- 

When an infinitesimal force is applied for positive times, one has 

d- 



vit) = —x(t) oc / C^{t')dt' . (4) 
dt Jo 

A straigthforward consequence of the above relations and of the simple identity 

{x\t))= f' f'c,{ti-t2)dtidt2 , (5) 
Jo Jo 

which suggests 

^ = f^dt oc ~ t^" , (6) 



in analogy with ([3]). On the other hand, it can be seen that such a formal argument 
is not rigorous and the actual scenario may become rather subtle, see e.g. [5]. For a 
detailed discussion the reader can see [1]. 
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In this paper we discuss the subdiffusive situation. Some works show that in such 
a case the expected resuh ([6]) seems to hold [6l [7]. This has been exphcitely proved 
in systems described by a fractional-Fokker-Planck [7j equation, where a generahzed 
Einstein relation has been shown {F is the perturbing force). 

lF(a;2(t)) 



x(t) = . (7) 

Models based on fractional Fokker-Planck equations, although interesting, usually 
are not directly derived from specific real systems; we therefore wondered whether a 
relation similar to ([7]) holds in more realistic models, such as in single file diffusion [8], 
which is a sub-diffusive system having many realizations in nature (e.g. transport 
in nanopores or narrow channels and zeolites, as well as car traffic on single lanes, 
pedestrian dynamics, etc.). The model used here consists of a one-dimensional gas 
of inelastic hard particles, moving on a large ring. To ensure a stationary state, 
particles exchange energy with an external thermostat. Tuning the characteristic 
time of the thermostat, the average volume fraction occupied by the gas and the 
restitution coefficient (from elastic to completely anelastic), one may observe a wide 
range of different stationary states, from a homogeneous density with Gaussian velocity 
distribution to strongly inhomogeneous spatial arrangement (clustering) with non- 
Gaussian statistics of velocities O |10l HH |T2[ |13] . Other authors have studied diffusion 
in granular gases without any external driving: in this case the gas is non-stationary 
(cooling regime) and one finds non-trivial exponents for diffusion [HI IT5j . 

The aim of this paper is to discuss the consequences of both sub diffusion and 
inelasticity in the more general context of linear response theory for statistically 
stationary states [HI |T7j. Let us briefly remind some general results [1]. Consider 
a dynamical system X(0) — >■ X()f:) = t/*X(0) whose time evolution can also be not 
completely deterministic {e.g. stochastic differential equations), with states X belonging 
to a A^-dimensional vector space. We assume a) the existence of an invariant probability 
distribution p(X), for which an "absolute continuity" condition is required (see jl] for 
details), and b) the mixing character of the system (from which its ergodicity follows). 
In our stochastic model the two above requests hold. Under these hypotheses, it is 
possible to derive (for details see [T6l[T71ll]) the following generahzed FD relation, valid 
when considering the perturbation at time of a coordinate Xj-. 



In the case of thermostatted Hamiltonian systems, on the other side, one has that 
p(q, p) oc exp(— /9?i(q, p)). From formula ([H]), therefore, one has that 

With a slight abuse of terminology, we will use the form "Einstein relation" to denote the 
time dependent Eq. ([9]). Let us note that its validity is a consequence of the Gaussian 
statistics of the velocity and the factorization of the stationary probability distribution. 
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i.e. positions and velocities are independent. In non-Hamiltonian systems, the shape of 
p(x) is not known in general, therefore ([8]) does not give a straightforward information. 
Nevertheless it can be exploited to get an interpretation of the results of a linear response 
experiment. We will analyze the response to small perturbations in the stationary state 
of a one- dimensional granular gas, discussing the response properties of the stationary 
state with its many "anomalies" with respect to an equilibrium state. 

We stress that the regimes considered here are always ergodic: this is a relevant 
difference with respect to the studies on the violations of the Fluctuation-Response 
relation, which considered glassy systems in the non-ergodic (aging) phase [18]. 

The model The model considered here consists of a gas of inelastic hard rods of 
mass 1, of linear size d, moving on a ring of length L. The rods interact also with a 
heating bath which mimics the effect of an irregular vibration injecting energy in the 
system. Until a collision occurs, the position Xi and the velocity Vi of i-th. rod obeys the 
following equations: 



where rjiit) is a Gaussian white noise with {rjiit)) = and {'i]i(t)'nj{t')) = Sij6{t — t'). 
When two rods i and j come into contact, their velocities Vi and Vj are instantaneously 
changed into v[ and f '• with the following rule: 



The meanings of r;, and are those of a typical thermalization time and a temperature, 
respectively, obtained if the system is elastic (r = 1). The coefficient of restitution 
r G [0, 1] determines the degree of inelasticity: after a collision, a fraction proportional 
to 1 — of the relative kinetic energy (i.e. kinetic energy in the center of mass frame) 
of the two particles is lost. When the particles are homogeneously distributed along the 
ring, the mean free path is given byA = 1/n — d = (1 — 0)/n where n = N/L is the 
number density and (f) = nd is the occupied volume fraction. The mean free time Tc 
is roughly estimated as X/^Jt^. In the rest of the paper we will tune n or r^, keeping 
fixed Tft = 1, in order to change the ratio between characteristic times a = Tc/ti,. For 
any value of r or a, the system reaches a statistically stationary regime where a kinetic 
temperature, denoted as "granular temperature", Tg = (t>^) can be measured. When 
a ^ 1, the coupling with the thermostat dominates the dynamics of the rods: they 
therefore remain thermalized and the system results at equilibrium at temperature Tf,: 
only spatial (rod-rod) correlations are expected at equilibrium, while velocities are not 
correlated, i.e. the global phase space probability distribution function (pdf) factorizes 
as: 




dt 





(10) 




(11) 



N 




(12) 



i=l 



FD in sub -diffusive systems 



5 



with Pv{v) a Gaussian distribution with variance T^. On the contrary, when a <^ 1, 
the effect of inelastic colhsion is strong enough to draw the system in a non-equihbrium 
stationary state whose properties are known from previous studies [9l [101 EH ESI fT3] . 
Non-Gaussian single particle velocity distributions and correlations among velocities and 
positions are the most relevant; these anomalies with respect to equilibrium become more 
and more pronounced as a or r are reduced. As a matter of fact, in this regime it is not 
correct to assume a factorization of the kind of Eq. f|T2l) . and the single particle velocity 
distribution, which is non-Gaussian, represents only a projection on a single degree 
of freedom of the full phase-space measure. We will see that the non-gaussianity of 
velocities is by far less important that the lack of factorization, which becomes relevant 
when the system is not dilute enough and which makes the Einstein relation ([9]) fail. 




Figure 1. Plot of the normalized autocorrelation C{t) versus time, for two cases, one 
elastic (full line) and the other inelastic (dashed line). In the left inset we show a 
blow-up of the exponential decay at early times. In the right inset you can find a blow 
up in log- log scale of the negative tail, together with a power law decay t~^^'^. Here 
= 0.1, and a « 0.9. 



The velocity autocorrelation function In Figure [T] we show the normalized 
autocorrelation function: C{t) = C^it) / 0^(0) = {v{t)v(fS)) /Tg for the velocity of a 
tagged particle (a tracer with the same properties of other particles). In both elastic and 
inelastic experiments, C{t) presents three main features: a) an exponential decay at early 
times, b) a negative minimum and c) asymptotically a power-law decay C{t) ~ —t^^/"^. 
The negative minimum is necessary to have subdiffusion, i.e. D = J^C{t) = 0, while 
the final power- law decay with 3/2 exponent is necessary to have (x^(t)) ~ t^^^. The 
initial exponential decay C{t) ~ exp( a more subtle nature. In ID 

one can argue that the tracer "discovers" the geometrical contraint after a long time. 
However, calculations based on collisions between non-correlated particles lead to wrong 
predictions for tcorr- Since this point is not closely related to the FD relation, we do not 
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discuss it in detail. Here we do not show the mean squared displacement as a function 
of time, already detailed in |19j: however the single-file diffusion scenario ~ t^/^ 

holds for any value of r, a and 0. 




Figure 2. Left: parametric plot of response R{t) versus normalized autocorrelation 
C{t). The dashed hne marks the Einstein relation R = C. Where not specified, 
Tb = 1. Right: ~C{t) and —R{t) versus t for elastic and inelastic cases at late times, 
with = 0.1 and Ti, — 1. 



The response to an impulsive perturbation The response to an impulsive perturbation 
is shown in Figure [2] for some choices of parameters. We have used a standard recipe 
to have a clean measure of response [20]: the system is let thermalize, then at time 
to is cloned. The original system evolves without perturbation, the copy is perturbed, 
i.e. the tagged tracer receives a small kick v{to) — > f'(to) = "^(^0) + with 6v <^ ^JTg 
to ensure linearity of the response. Then the copy is evolved using the same noise 
realization as for the original system and the response is given by the dynamical average 
R{t) = {v'(to + t) — v{to + t))/6v over many realizations of the experiment. In Figure [2] 
we show representative cases where the Einstein relation R{t) = C{t) is verified within 
numerical precision. This happens for elastic cases, or cases at low inelasticity 1 — r ^ 1 
and low packing fraction, and also for cases at high inelasticity, provided that u ^ Tc. 
This last setup corresponds to a very fast action of the thermal bath which practically 
removes the effects of inelastic collisions. Similar results have been obtained, previously, 
for 2d driven granular gases [211 ESI EHl EH [25] . As shown in the top-right frame, for the 
elastic case, the relation R{t) = C (t) is fairly verified also at late times in the power- law 
tail. The inelastic case (see bottom-right frame) displays a small violation at such large 
times: note that this small violation corresponds to both R{t) and C{t) very close to 
zero. 

In Figure [3] the parametric plot of response versus correlation is displayed for 
cases where the Einstein relation is no more verified. The departure from the equality 
R(t) = C (t) can be quite strong: it increases with the packing fraction 0, the inelasticity 
1 — r and the rescaled bath time Th/Tc = 1/a. In all cases we observe R{t) < C{t). In 
Fig. [3] we have stressed the dependence on a, which can be tuned changing Tb at fixed 
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Figure 3. Parametric plot of response R{t) versus normalized autocorrelation C{t). 
The dashed line is the Einstein relation R = C. All data are obtained "with restitution 
coefficient r — 0.6. On the left: the packing fraction is constant (j) — 0-1 and is 
changed, resulting in different values of Tc- The ratio a = Tc/t;, is given for simplicity. 
On the right: Tf, = 1 is kept constant, "while <j) is changed. In the insets the correlator 
C„2 „2, discussed in the text, is displayed as a function of the varying parameter, for 
elastic and inelastic systems. 



r and (p. In all experiments we have verified to be in the linear response regime. 

Origin of the violation of the Einstein relation As anticipated in the description of the 
model, and in agreement with the observation done in [25], the Einstein relation no more 
holds when the factorization of the phase-space pdf expressed by Eq. ( fT2l) is violated. 
For reasons of space we do not show the probability density function of one-particle 
velocities, which are not far from the Maxwell-Boltzmann distribution. Violations of 
Gaussianity have been shown in [25] to be not relevant for the FD relation, because 
autocorrelations at different orders are almost proportional, i.e. {v{0)v{t)) / {v"^) ^ 
{v{Oyv{t))/{\v\^) ^ (f (0)^f etc. This is confirmed by Direct Monte Carlo 
simulations, where an almost perfect factorization of the degrees of freedom in the 
phase-space pdf is satisfied: in such simulations, even with a stronger departure from 
Gaussianity, the Einstein relation always holds. 

Many ways of characterizing the breakdown of phase-space factorization can be 
employed. A simple one is displayed in the inset of Fig. [3j 

where Svf = vf — Tg. When 0^2 > 0, the squared velocities of two adjacent particles 
are correlated. It is evident that this correlation increases when a is decreased. The 
same is observed tuning the other parameters, such as decreasing r or increasing (p. 



Conclusions Drawing the conclusions, we stress the twofold nature of this study. On 
one side, for the elastic single-file diffusion, which is a less abstract model than fractional 
Fokker-Planck, we have obtained a good agreement between R{t) and C{t), in all time 
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ranges, confirming the validity of the FD ("Einstein") relation. On the other side 
we have explored the effects of inelasticity: in this case one has a non-equilibrium 
stationary state where strong correlations among different particles are present, therefore 
the factorization (fT2|) fails and only a more general FD relation ([8]) holds. At small 
inelasticity, small packing fraction and/or for fast thermostats, the Einstein relation 
is recovered, because the lack of factorization is weak, as previously observed in 2d 
granular gases [211 l22l ISHl l23l l2ll [25]. A quantitative characterization of the departure 
from factorization is under investigation, with the aim of proposing, as a first step, a 
joint two-particles (first neighbours) velocity distribution: we expect to obtain, from 
this study, a first explicit correction formula to the Einstein relation. 
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